Phonon-assisted optical absorption in silicon from first principles 
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The phonon-assisted interband optical absorption spectrum of silicon is calculated at the quasi- 
particle level entirely from first principles. We make use of the Wannier interpolation formalism 
to determine the quasiparticle energies, as well as the optical transition and electron-phonon cou- 
pling matrix elements, on fine grids in the Brillouin zone. The calculated spectrum near the onset 
of indirect absorption is in very good agreement with experimental measurements for a range of 
temperatures. Moreover, our method can accurately determine the optical absorption spectrum of 
silicon in the visible range, an important process for optoelectronic and photovoltaic applications 
that cannot be addressed with simple models. The computational formalism is quite general and 
can be used to understand the phonon-assisted absorption processes in general. 



The phonon-assisted absorption of light in materials is 
an important optical process both from a fundamental 
and from a technological point of view. Intraband light 
absorption by free carriers in metals and doped semicon- 
ductors requires the additional momentum provided by 
the lattice vibrations or defects, while phonon-assisted 
processes determine the onset of absorption in indirect- 
band-gap semiconductors (Fig. [I]). Moreover, the value 
of the direct band gap in silicon (3.4 eV[JJ) is large and 
precludes optical absorption in the visible. However, sil- 
icon is a commercially successful photovoltaic material 
because of the indirect optical transitions that enable 
photon capture in the spectral region between the in- 
direct (1.1 eV[2]) and direct band gaps. 

Despite their importance, at present, only a very lim- 
ited number of first-principles studies of phonon-assisted 
optical absorption spectra exist. Ab initio calculations of 
direct optical absorption spectra including excitonic ef- 
fects have already been performed for Si and other bulk 
semiconductors [Jj and the underlying methodology is 
presently well established [H [5] . Phonon-assisted absorp- 
tion studies are more involved, however, and the associ- 
ated computational cost is much higher than the direct 
case. The calculation of the indirect absorption coef- 
ficient involves a double sum over fc-points in the first 
Brillouin zone (BZ) to account for all initial and final elec- 
tron states. In addition, these sums must be performed 
with a very fine sampling of the zone to get an adequate 
spectral resolution. The computational cost associated 
with these BZ sums is in fact prohibitive with the usual 
methods. Phonon-assisted absorption calculations have 
been done for the special case of free-carrier absorption 
in semiconductors [B], where the carriers are initially lim- 
ited to a region near the Y point of the first BZ, but a 
full calculation using brute-force methods for the general 
case remains beyond the reach of modern computers. 

The difficulty of zone-integral convergence can be 



addressed with the maximally-localized Wannier func- 
tion interpolation mcthod[7]. Using this technique, the 
quasiparticle energies [8] and optical transition matrix 
elements [S] can be interpolated to arbitrary points in the 
BZ at a minimal computational cost. Moreover, this in- 
terpolation method has been generalized[10 to obtain 
the electron-phonon coupling matrix elements between 
arbitrary pairs of points in the first BZ. 

In this Letter, we developed a first-principles computa- 
tional method, based on a Wannier-Fourier interpolation 
formalism, to calculate the phonon-assisted optical ab- 
sorption spectra of materials from first principles and ap- 
plied it to the case of interband absorption in silicon. The 
calculated spectra near the absorption onset are in very 
good agreement with experimental results for a range of 
temperatures. Moreover, we were able to reproduce the 
absorption spectrum in the energy range between the in- 
direct and direct band gaps (1.1 - 3.4 eV), a spectral re- 
gion that cannot be accessed by standard model calcula- 
tions. This region covers the entire visible spectrum and 
is important for optoelectronic applications. The com- 
putational formalism is quite general and can be used to 
predict and analyze the phonon-assisted optical absorp- 
tion spectrum of any material. 

To calculate the phonon-assisted absorption coeffi- 
cient, we use the Fermi's golden rule expression [6j ITT] : 

a(w) =2-^^-L 1 J_- £ \\-{S l + S 2 )\ 2 
wcn r u Veen N k N q 

X Pd(e jt k+ q — e,fe - fojj ± ftw V q), (1) 

where fku and A are the energy and polarization of the 
photon and n r (u>) is the refractive index of the material 
at frequency uj. The generalized optical matrix elements, 
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Si and S2, are given by 
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and correspond to the two possible paths of the indi- 
rect absorption process (Fig. [T]). They are determined in 
terms of the velocity (v) and electron-phonon coupling 
(g) matrix elements, as well as the real (e n k) and imag- 
inary (r n fc) parts of the quasiparticle self-energies. The 
factor P accounts for the carrier and phonon statistics, 
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The upper (lower) sign corresponds to phonon emission 
(absorption). 

The Kohn-Sham eigenvalues were calculated within 
the local density approximation (LDA)[12] to density 
functional theory using a plane-wave pseudopotential 
approach [T3] with a kinetic energy cutoff of 35 Ry. The 
ground state charge density was determined on a BZ grid 
of 14x14x14 fe-points. Quasiparticle energies within the 
GW approximation for the self-energy operator [14] were 
determined on a 6x6x6 grid and interpolated through- 
out the BZ through the use of the maximally-localized 
Wannier function formalism [7 . We included 34 elec- 
tronic bands in the coarse-grid calculation and extracted 
26 Wannier functions, which reproduce the LDA band- 
structure 10 eV below and 30 eV above the Fermi level. 
The interpolated quasiparticle band structure of silicon is 
shown in Figure [I] The indirect (1.3 eV) and direct (3.3 
eV) quasiparticle band gaps are in good agreement with 
previous calculations [2] and experiment. The same for- 
malism has been used to interpolate the velocity matrix 
element s [31 H], including the renormalization required [15] 
after the GW corrections. The real (ei) and imaginary 
(£2) parts of the dielectric function and the refractive 
index due to direct transitions, required in Eq. [I] to de- 
termine the absorption coefficient, were also determined 
at the quasiparticle level for a range of photon frequen- 
cies. Lattice dynamics are calculated using density func- 
tional perturbation theory [IB"]. The electron-phonon cou- 
pling matrix elements are calculated on the same coarse 
grid of electronic points, while the dynamical matrices 
and phonon-potential perturbations are calculated on a 
6x6x6 grid of momentum-space vectors [TU] and interpo- 
lated for arbitrary pairs of points in the first BZ using 
the EPW codejTT]. For the calculations of the velocity 
and electron-phonon coupling matrix elements we used 
the LDA wave functions, because their overlap with the 
GW-corrected ones is better than 99.9% [HJ. 

Phonon-assisted optical absorption in indirect-band- 
gap semiconductors occurs for photons with energies 
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FIG. 1. (a) Quasiparticle band structure of silicon calcu- 
lated within the GW approximation and interpolated with the 
Wannier formalism. The arrows indicate the lowest-energy 
phonon-assisted optical absorption processes across the indi- 
rect band gap. Solid lines denote optical transitions, while 
dashed lines correspond to electron-phonon scattering events. 
The two terms, Si and S2, of Eq.[l]are represented by paths 
(1) and (2) respectively, (b) Density of electronic states versus 
the quasiparticle energy. The density of states of the occupied 
bands has been highlighted. 



greater than the indirect band gap minus (plus) the en- 
ergy of the phonon absorbed (emitted) to assist the tran- 
sition. The onset of indirect absorption is calculated over 
a wide range of temperatures in bulk silicon through Eq.[l] 
and the results are shown in Fig. [2] Each curve displays 
a characteristic knee, arising from the different energy 
onsets of the phonon-absorption and phonon-emission 
terms, which becomes smoother with increasing temper- 
ature. The calculated data are in good agreement with 
experimental results 18J for all temperatures measured. 
For these calculations, we used fine grids of 40x40x40 
for the k and q sums in Eq. [T] respectively. These fine 
grids yield converged optical spectra with an energy res- 
olution of 14 meV, which is quite small and necessary to 
resolve the fine features near the absorption onset. Al- 
though the experimental data near the edge can be fit 
with simple parameterized forms [T5], to our knowledge 
they have not been calculated entirely from first princi- 
ples previously. 

In addition to the absorption onset, we are interested 
in the phonon-assisted absorption spectrum in the en- 
ergy range between the indirect and direct band gaps, 
covering the visible range. This spectral region involves 
transitions between valence and conduction band states 
away from the band extrema and, as a consequence, can- 
not be modeled with simple parameterized forms. On the 
contrary, because of the large number of electronic states 
and phonon modes involved, first-principles calculations 
are the only computational tool that can access this spec- 
tral region. The interpolation of the ab initio quantities 
within the Wannier-Fourier formalism makes the calcula- 
tion feasible on modern computers. The calculated spec- 
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FIG. 2. Onset of the phonon-assisted optical absorption in silicon, as a function of photon energy and temperature. The 
theoretical results (solid lines) are in good agreement with experiment (dashed lines) . Experimental data are from Ref. 1181 The 
theoretical curves have been shifted horizontally to match the onset of the experimental spectra. 
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FIG. 3. Calculated (solid lines) and experimental (circles) 
absorption coefficient of silicon in the energy range between 
the indirect and direct gaps, for two temperatures. Experi- 
mental data are from Ref. I19I The theoretical spectra have 
been shifted to match the experimental absorption onset. 



tra with an energy resolution of 30 meV (Fig. [3]) converge 
with zone-sums of 24x24x24 electronic and 24x24x24 
phonon points. The imaginary part of the electron self- 
energy for the intermediate electronic states was set to a 
constant value (100 meV). However, the calculated data 
are not very sensitive to the particular value of the quasi- 
particle lifetime for photon energies in this spectral re- 
gion. 

To facilitate comparison with experiment, the theoret- 
ical absorption spectra of Figs. [2] and [3] have been rigidly 
shifted to the left along the energy axis by 0.15-0.23 eV in 
order to match the onset of the experimental absorption 
curves. This shift is needed to account for the numerical 
difference between the calculated and experimental band 
gap, and for finite-temperature effects on the quasipar- 
ticle energies which we have not considered explicitly. 



Although the GW method is the most accurate first- 
principles computational formalism for the calculation of 
quasiparticle properties presently available, it typically 
yields absolute quasiparticle energies accurate to 0.1 eV. 
In our particular case, we found that the calculated band 
gap is also within this error bar larger than the exper- 
imentally measured value. We note that no other first- 
principles method is presently available to guarantee a 
more accurate result. Moreover, the only temperature 
dependence we considered in our calculations is for the 
phonon occupation numbers. However, the quasiparticle 
energies themselves are temperature-dependent because 
of additional finite-temperature effects, such as lattice 
expansion and electron-phonon renormalization [TU] |2"0T - 
[2^] . We found that the thermal-expansion correction to 
the indirect band gap is small and amounts only to a 2.5 
meV increase of the LDA band gap as the lattice con- 
stant increases from the K to the 400 K value [25], in 
agreement with Ref. On the other hand, empirical 
pseudopotential calculations have shown that electron- 
phonon renormalization effects are stronger and decrease 
the band gap of silicon by approximately 0.05-0.1 eV for 
temperatures in the range 0-400 K [22] . We note that the 
determination of electron-phonon-coupling corrections to 
quasiparticle energies from first principles is still a sub- 
ject of ongoing research [53]. The cumulative effects of 
this electron-phonon band-gap renormalization (0.05-0.1 
eV) and the intrinsic accuracy of the GW method (or- 
der of 0.1 eV) explain the difference between the onsets 
of the theoretical and experimental data that we need to 
take into account when comparing our calculated spectra 
to experiment. Moreover, this electron-phonon coupling 
correction to the quasiparticle energies may have an ef- 
fect on the shape of the absorption spectra near the onset 
of indirect transitions in Fig. [2] 

Although excitonic effects, arising from the electron- 



4 



hole interaction, are potentially important for optical 
processes and in general affect the direct absorption spec- 
tra even for photon energies far from the absorption 
edge[l], we found that they are not as important for the 
case of indirect absorption. Sharp features that appear 
near the onset of indirect absorption are attributed to 
excitonic effects [TH1 126) . Our calculations, however, are 
based on quasiparticle theory and do not account for the 
electron-hole interaction that gives rise to these excitonic 
features. Nevertheless, the calculated absorption spectra 
are in very good quantitative agreement with experimen- 
tal data, pointing to a weaker role of the electron-hole 
Coulomb interaction for the case of indirect optical tran- 
sitions. This is probably because the band-extrema wave- 
functions in indirect-gap materials are located at differ- 
ent points of the BZ and hence the wavefunction overlap, 
which determines the magnitude of the Coulomb inter- 
action between them, is small. Therefore, the phonon- 
assisted spectra can to a large extend be explained at 
the quasiparticle level of theory, without the need to ac- 
count for excitonic effects. 

The computational formalism we developed is based 
on first-principles methods and can be used to study 
the fundamental physics of phonon-assisted absorption 
in materials in general. It can complement experimental 
studies to shed light on the microscopic phonon-assisted 
transition mechanisms and address questions that are 
not accessible by experimental techniques. Moreover, the 
method can be used to analyze the phonon-mediated op- 
tical properties of technologically important materials for 
optoelectronic applications. E.g., it can investigate the 
role of phonon-assisted optical processes in silicon pho- 
tonics, or it can predict the photovoltaic performance of 
indircct-band-gap materials. 

In conclusion, we used a Wannier-Fourier interpola- 
tion technique to calculate the phonon-assisted optical 
absorption spectra of silicon at the quasiparticle level. 
The calculated spectra are in very good agreement with 
experimental measurements, both near the absorption 
onset and in the spectral region between the indirect and 
direct band gaps for any lattice temperature. The first- 
principles computational formalism is very general and 
can be used to study the fundamental physics of phonon- 
assisted absorption, as well as the phonon-mediated op- 
tical properties of optoelectronic materials. 
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